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FINAL  REPORT 


ARO  PROPOSAL  NUMBER:  21078-GS 

TITLE  OF  PROPOSAL:  Fundamental  Studies  on  Basin  Scale  Hydrology 
CONTRACT  NUMBER:  DAAG29-84-K-0064 

NAME  OF  INSTITUTIONS:  University  of  Mississippi  and 

Oregon  State  University 

AUTHORS  OF  REPORT:  V.  K.  Gupta  and  Ed  Waymire 

1.  OBJECTIVES 


The  proposed  research  had  two  broad  objectives.  These  objectives,  as 
stated  in  the  proposal,  were: 

(I)  Identify  the  ' geomorphologic  laws'  of  drainage  networks  in  three 
dimensions  (i.e.,  taking  into  account  both  elevations  and 
horizontal  spatial  dimensions)  which  determine  the  path  of  runoff 
to  the  basin  outlet  fed  into  the  channel  network  from  the  adjacent 
hillslope  regions  and  (b)  undertake  a  theoretical  representation 
of  these  'laws'  on  the  basis  of  a  few  general  postulates. 

(II)  A  coupling  of  the  stochastic  rainfall  field  in  space-time  with  the 
basin  scale  response  for  investigating  the  statistical  behavior  of 
the  streamflow  process. 

Dr.  Gupta  was  the  Principal  Investigator  in  charge  of  the  research  on 
Objective  I,  and  Dr.  Waymire,  as  a  Sub-Contractor,  on  Objective  II. 
However,  both  Drs.  Gupta  and  Waymire  worked  very  closely  on  the  entire 
research  effort. 


2.  MAIN  FINDINGS  AND  RESULTS 

This  entire  research  effort  was  built  on  two  mathematical  themes 
which  are  necessary  for  a  fundamental  theoretical  understanding  of 
climate-hydrology-geomorphology  interactions  and  for  solving  the 
practical  problem  of  predictions  from  ungaged  basins.  These  themes 
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include  (i)  an  ensemble  view  or  a  probabilistic  framework,  and  (ii) 
self-similarity  and  scaling  in  space  and/or  time.  The  following  pages 
briefly  describe  the  main  findings  and  progress  separately  under  each 
objective.  However,  throughout  this  research,  our  directions  and 
progress  on  either  one  of  the  two  objectives  greatly  influenced 
directions  and  progress  on  the  other. 
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magnitude  m  specifies  the  systematic  variability.  Theoretical 
predictions  were  tested  against  empirical  observations  by  predicting 
a  network  structure  function  called  the  link  concentration  function 
(lcf).  The  lcf,  denoting  the  number  of  links  as  a  function  of 
elevation,  was  introduced  by  Gupta,  Waymire  and  Rodriguez- I turbe 
[1986].  Four  alternate  sets  of  assumptions  on  p(m)  were  tested 
using  data  from  four  basins  from  different  climates  and  geologies. 
Rather  unexpectedly,  these  tests  led  us  to  the  discovery  that  link 
height  dfs  exhibit  invariance  under  magnification  or  reduction  of  the 
scale  parameter .  This  invariance  property  is  referred  to  as 
'statistical  self-similarity’.  It  provides  the  following 

mathematical  representation  of  the  random  height  of  a  link  of 
magnitude  m  in  a  channel  network: 

H(m)  i  m~S  H(l)  ,  (1) 

where  6  is  a  fundamental  scaling  exponent,  and  the  symbol  '  =  ' 
in  eq.  1  means  that  the  random  variables  on  both  sides  of  it  have  the 
same  dfs.  Self-similarity  provides  a  fundamental  theoretical 

explanation  of  some  purely  empirical  relationships  regarding  mean 
gradients  and  other  geometries  in  channel  networks.  Preliminary 
formulations  and  tests  of  the  above  findings  are  given  in  Mesa  [1986] 
and  Gupta  and  Mesa  [1988],  Detailed  theory  that  led  to  this 
discovery  is  contained  in  Gupta  and  Waymire  [1988], 
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2.1.1.  Generalizations  of  Self-Similarity  and  Prediction  From  Uneaged 
Basins 

We  will  briefly  describe  the  issues  concerning  theoretical 

generalizations  of  self-similarity  in  link  heights  and  gradients  for 

channel  network  structures  and  peak  flows  in  river  basins.  These 

generalizations  have  a  direct  bearing  on  predictions  from  ungaged 

basins.  More  details  are  given  in  Gupta  and  Waymire  [1988], 

Under  the  physically  reasonable  assumption  of  statistical 

independence  between  link  gradients  and  link  lengths,  and  the  random 

model  postulate  on  link  lengths,  it  follows  that  link  gradients  in  a 

channel  network  are  also  self-similar.  The  scaling  exponents  8 

for  link  heights  and  link  gradients  are  the  same.  Moreover,  this 

exponent  is  unique .  Let  s^m)  denote  the  a-quantile  of  a  link 

gradient  S(m)  of  magnitude  m,  i.e.,  P(S(m)  >  s  (m) )  =  a.  It 

a 

follows  from  self-similarity  that, 

s  (m)  =  s  (1)  m~B  (2) 

oc  a 

Eq.  2  leads  to  a  fundamental  theoretical  explanation  that,  purely 
empirical  but  widely  observed,  scaling  of  the  means  of  link 
gradients,  widths,  etc.,  reported  in  the  literature  is  due  to  self¬ 
similarity  in  these  variables. 

Beyond  the  means,  the  empirical  relationship  for  quantiles  of 
channel  gradients  and  the  theoretical  relationship  (2)  differ  in  an 
important  way.  The  empirically  computed  exponents  for  gradients  are 
observed  to  vary  with  the  frequency  a,  whereas  the  theoretical 
exponent  8  in  (2)  cannot  exhibit  such  a  variation.  Therefore, 
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tests  of  self-similarity  in  channel  gradients  and  other  geometries 
and  peak  flows  were  explicitly  carried  out  in  an  M.S.  Thesis  by 
Cadavid  [1988],  These  tests  clearly  show  that  self-similarity  is 
inadequate  to  describe  the  spatial  variability  in  channel  geometries 
and  peak  flows  in  a  channel  network.  This  feature  would  come  about 
if  link  gradients  in  a  network  are  not  shaped  by  a  single  dominant 
flow  frequency,  but  rather  a  whole  spectrum  of  frequencies. 

A  theoretical  generalization  of  self-similarity  as  a  mul ti- 
scaling  invariance  is  needed  before  it  can  be  tested  from  data  for 
peak  flows  and  link  gradients  in  channel  networks .  A  preliminary 
analysis  of  this  kind  for  peak  flows  was  carried  out  in  an  M.S. 
thesis  by  Cadavid  [1988]  with  very  encouraging  results.  However, 
the  mathematical  foundations  of  the  multiscaling  construct  used  in 
that  thesis  remain  quite  obscure.  A  precise  formulation  of  multi¬ 
scaling  is  not  known  and  will  require  a  lot  more  work.  Should  this 
line  of  investigation  prove  successful  in  the  future,  it  would 
provide  exciting  new  possibilities  unthought  of  before  towards  flood 
predictions  from  ungaged  basins  and  regionalization  of  flood 
frequency  analysis. 
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Therefore,  an  understanding  of  such  theoretical  constructs  becomes 
mandatory  if  we  are  to  be  able  to  understand  how  rainfall  fields 
direct  space-time  variability  in  streamflows.  Keeping  in  view  this 
focus,  we  will  only  summarize  those  results  on  rainfall  which  dealt 
with  the  scaling/self-similarity  issues. 

Self-similarity  can  arise  most  naturally  as  an  asymptotic 
property  under  a  wide  variety  of  mathematical  constructs.  The  ideas 
on  scaling  pertinent  to  rainfall  were  explained  and  illustrated  by 
Waymire  [1985]  in  an  expository  article. 

Analysis  of  spatial  rainfall  from  GATE  data  sets  shows  a  scaling 

property  in  the  second  moment  of  spatially  averaged  rainfall.  The 

2 

averages  are  taken  over  regions  ranging  from  4  x  4  km  to  200  x  200 
2 

km  .  The  same  data  sets  have  also  been  used  to  model  the  probability 
density  of  spatially  averaged  rainfall  by  a  log-normal  density.  In 
Waymire  and  Gupta  [1988]  it  was  shown  that  self-similarity,  defined 
in  the  usual  sense  of  distribution  functions ,  and  lognormality  are 
mutually  inconsistent.  This  suggests  modifications  of  self¬ 
similarity  or  lognormality,  or  both,  for  modeling  spatial  rainfall. 
Indeed,  multi-scaling  ideas  are  being  explored  by  us  and  others  in 
describingh  space-time  rainfall. 

Finally,  in  an  article  by  Gupta  and  Waymire  [1987]  it  was  shown 
that  if  a  rainfall  field  is  self-similar  in  space,  then  this  property 
carries  over  in  time  via  Taylor's  hypothesis.  This  means  that 
temporal  rainfall  becomes  self-similar.  However,  the  presence  of 
dissipation  in  rainfall  destroys  Taylor's  hypothesis  after  a  time 
span  of  the  order  of  a  typical  life  cycle  of  a  convective  rainfall 
cell.  This,  in  turn,  leads  to  the  conclusion  that  self-similarity  in 
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temporal  rainfall  can  manifest  at  most  over  a  short  time  span  but  not 
beyond  that. 
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